Collective Properties of X-ray Binary Populations of Galaxies III. 
The Low-mass X-ray Binary Luminosity Function 

H. Bhadkamkar^'^ and P. Ghosh^ 
Received ; accepted 



^Department of Astronomy & Astrophysics, Tata Institute of Fundamental Research, 
Mumbai 400 005, India 

^Astronomy & Astrophysics, Raman Research Institute, Bengaluru, 560080 



- 2 - 
ABSTRACT 

Continuing our exploration of the collective properties of low-mass X-ray bi- 
naries (LMXBs) in the stellar fields (i.e., outside globular clusters) of normal 
galaxies, we compute in this paper (Paper III in the series) the expected X-ray 
luminosity function (XLF) of LMXBs, starting from the results of the previous 
paper in the series (Paper 11). For this, we follow the evolution of the pre-LMXB 
population (described in Paper II) through Roche-lobe contact and the conse- 
quent LMXB turn-on, LMXB evolution through the accretion phase, and the 
eventual conclusion of accretion and LMXB turn-off. We treat separately two 
classes of LMXB evolution, the first (a) being close systems whose initial orbital 
periods are below the bifurcation period, wherein the companion is on the main 
sequence when Roche-lobe contact occurs, the subsequent evolution is driven 
by angular-momentum loss from the system, and the second (b) being wider 
systems whose initial orbital periods are above the bifurcation period, wherein 
the companion is on the giant branch when Roche-lobe contact occurs, and the 
evolution of these systems is driven by the nuclear evolution of the companion. 
We obtain model luminosity profiles L{t) for individual LMXBs of both classes, 
showing that they are in general agreement with those in previous works in the 
subject. We point out that the basic features of the luminosity profile in the 
angular-momentum-loss driven case can be well-understood from scaling laws 
first pointed out by Patterson (1984) in a related context, which we call Patter- 
son scaling. We then compute the LMXB XLF by "folding in" the inputs for 
the pre-LMXB collective properties and formation rates from Paper II with the 
above luminosity profiles. Because of the long timescale (~ 10^ yr) on which 
LMXBs evolve, this computation becomes more involved, as one needs to keep 
track of the evolution of the star-formation rate (SFR) on the same timescale. 
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and we use star-formation histories (SFH) given by canonical models. We com- 
pare the observed LMXB XLF with our computed one, keeping in mind that we 
have included only neutron-star systems in this work, so that there would be an 
unaccounted-for population of black-hole binaries at the high-luminosity end. We 
show that a qualitative similarity already exists between the two, and we discuss 
the role of the giant fraction, i.e., that fraction of all LMXBs which harbors a 
low-mass giant companion, on the shape of the XLF at the high-luminosity end. 

Subject headings: binaries: close - stars: evolution - stars: neutron - stars: low-mass 
- supernovae: general - X-rays: binaries - X-rays: galaxies 
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1. Introduction 

This is Paper III in our series of papers exploring the basic physics underlying the 
collective properties of accretion powered X-ray binaries. In Paper II, we described our 
study of the constraints on the formation of the so-called pre-low-mass X-ray binaries 
(pre-LMXBs) , which is a relatively rare process, and our calculation of the collective 
properties (i.e., distributions of their essential parameters) and the rates of formation of 
these pre-LMXBs. In this paper, we use the results of Paper II to compute the LMXB 
X-ray luminosity function (XLF). 

Our procedure here is to first follow the evolution of the pre-LMXB population through 
Roche-lobe contact and consequent LMXB turn-on, LMXB evolution through the accretion 
phase, and eventual conclusion of accretion and LMXB turn-off. To this end, we calculate 
the X-ray luminosity evolution of individual model systems in our computational scheme, 
and so construct luminosity profiles, L{t), for these systems. We treat separately the two 
well-known categories of LMXB evolution. The first is that for close systems whose initial 
orbital periods are below the so-called bifurcation period, wherein the companion is on the 
main sequence when Roche-lobe contact occurs and mass transfer starts. Evolution of 
such systems is driven by angular-momentum loss from the systems through gravitational 
radiation and magnetic braking, and their orbits shrink as mass transfer proceeds. The 
second is that for wider systems whose initial orbital periods are above the bifurcation 
period, wherein the companion has finished its main-sequence evolution and is on the 
giant branch when Roche-lobe contact occurs and mass transfer starts. Evolution of these 
systems is driven by nuclear evolution and expansion of the companion, and their orbits 
expand as mass transfer proceeds. We obtain model luminosity profiles for both categories, 
and show that they are in general agreement with those in previous works in the subject. 
We point out that the basic features of the luminosity profile in the angular-momentum-loss 



- 5 - 



driven category can be accounted for by scaling laws first pointed out by iPattersonI ( ll984j ) 
in a related context, which we call Patterson scaling. 

We obtain the LMXB XLF by combining the above inputs. Because of the long 
timecasle (~ 10^ yr) on which LMXBs evolve, the computation here becomes more involved, 
as one needs to keep track of the evolution of the star-formation rate (SFR) on the same 
timescale, since the rate of formation of primordial binaries, which evolve into LMXBs, is 
determined by the SFR. Accordingly, we construct a scheme for "folding in" the inputs from 
the pre-LMXB collective properties and formation rates, together with the star-formation 
history (SFH), i.e., the time- variation profile of the SFR, as given by canonical models 
in wide use. We thus obtain the theoretical LMXB XLF and study its properties and 
their variations with some parameters characterising the primordial binary distribution the 
supernova process. 

We compare the observed LMXB XLF with our computed one, keeping in mind 
that we have included only neutron-star systems in this work, so that there would be an 
unaccounted-for population of black-hole binaries at the high-luminosity end. We show that 
a qualitative similarity already exists between the two, and we discuss the role of the giant 
fraction, i.e., the fraction of all LMXB systems harboring a low-mass giant companion, on 
the shape of the XLF at the high-luminosity end. 

The rest of the paper is organized as follows. In Sec|2], we describe the essential aspects 
of LMXB evolution that we need for our purposes here. We give the prescriptions for 
angular momnetum loss by graviational radiation and magnetic braking, and we describe 
the two categories of LMXB evolution in terms of the bifurcation period. We obtain the 
luminosity profile for angular-momentum-loss driven evolution, and introduce and discuss 
the Patterson scaling (1984) which accounts for it. We then obtain the luminosity profile 
for evolution driven by nuclear evolution of the companion. In Secj3|, we describe our 
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scheme for computing the XLF through the device of convolving the inputs on pre-LMXBs 
from Paper II with these luminosity profiles, together with the SFH for taking into account 
the simultaneous evolution of the SFH. We discuss the properties of our computed XLF, 
exploring the effects of some primordial-binary parameters and supernova characteristics 
on it. In SecJH we compare our computed XLF with observations, and discuss the results. 
Our conclusions are presented in Sec J4.2[ 



2. LMXB Evolution 

A pre-LMXB system is initially detached. Mass transfer and accretion onto the neutron 
star is not possible in this state and becomes so only when the system evolves, its orbit 
shrinks due to loss of angular momentum, and Roche-lobe contact is achieved. The orbital 
angular momentum of the binary is given by 

J' = G^^a (1) 

Here Mc and M^s are the masses of the companion and the neutron star respectively 
and Mt = Mjv5 + is the total mass of the system, a is the binary separation. 

We describe binary orbital evolution due to angular momentum loss by differentiating 
Eq.(IT]) with respect to time. We first note that Mt = M^s + ^c- Since the companion is 
the mass donor, we can write Mc = —M where M is the mass-transfer rate. Allowing for 
non-conservative mass transfer in general, we write M^s = jM, i.e., a fraction 7 of the 
mass transferred by the donor is accreted onto the neutron star and the rest of the mass 
is lost from the system, reducing its total mass at a rate Mt = ~(1 ~ 7)M. This mass 
loss implies an appropriate loss of angular momentum from the system. It is customary to 
represent this loss of angular momentum in terms of the specific angular momentum of the 
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system (iBelczynski et al.ll2008l : iPfahl et al.ll2003l ). expressing the result as J / J = fiMt/Mt. 
Here, fi is typically of order unity, its exact value depending on the binary parameters 
as well as the evolutionary state of the donor. By inserting the expressions for the rates 
of change of various quantities as described above, the differential form of eqn. [1] can be 
rewritten as: 



J 



M 



NML 



7g - 1 + 



(l-7)(l + 2/i) q 
2 l + q 



+ 



2a 



(2) 



Here, {J/ J) nail is the rate of loss of angular momentum from the system due to 
processes other than the above mass loss from the system, and q = M^/Mns is the mass 
ratio. We now consider such processes of angular momentum loss. 



2.1. Mechanisms of angular momentum loss 

Two mechanisms have been widely studied in the literature that cause the angular 
momentum loss in the pre-LMXB systems. These are gravitational radiation and magnetic 
braking. We now briefly review these two mechanisms and discuss their effects on the 
binary parameters. 

It is well-known that two point masses orbiting around ea ch other would emi t 



gravitationa l waves according to the general theory of relativity (jPeters fc Mathews 



1963 



Peters 



196J), and the rate of evolution of the orbital separation (as also the eccentricity 



i n cas e of eccentric orbits) due to this effect is well-known. It was shown by 



Faulkner 



(jl97ll ) that this effect can be important in close binary systems like CVs and LMXBs. 
The timescale for the rate of angular momentum loss is comparable to typical LMXB 
evolutionary timescales, which makes it an essential effect in LMXB evolution. The rate of 
loss of the angular momentum due to this effect is given for a circular orbit by: 
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- j = -0.831 ' * (3) 

- / GW ^ 

We shall use the following system of units throughout this paper, unless stated otherwise. 
All the masses and the radii/ distances are given in terms of the solar mass and the solar 
radius respectively, and time is in the units of Gyr. X-ray luminosity of the LMXBs will he 
given in units o/ 10^^ erg/s and will he denotes as L^q. 

We observe the strong dependence on the orbital separation in Eq.Q, which 
makes gravitational radiation completely ineffective at large separations. However, since 
pre-LMXB orbits are already very compact due to large contractions during the CE phase, 
this mechanism can be, and often is, the dominant one determining the orbit shrinkage that 
converts pre-LMXBs into LMXBs. 

M agnetic braking wa s sugg ested as a possible mechanism for angular momentum 



loss by IVerbunt &: Zwaad (Il98ll ) in order to explain the observed high rates of the mass 
transfer in some LMXB systems. This mechanism depends upon the mass loss due to a 
magnetic wind from a tidally locked companion, working as follows. The magnetic field 
of the companion makes the matter in its coronal region co-rotate with the star, so that 
this matter has a large specific angular momentum, and consequently even a moderate 
mass-loss rate can lead to a large rate of angular momentum loss. Since the companion is 
tidally locked in pre-LMXBs and LMXBs, this loss of angular momentum will ultimately 
come from the orbital angular momentum of the binary, making it shrink. The rate of loss 
of angular momentum due to magnetic braking is given by 

A =-57.2„#^ (4) 

Here Rc is the radius of the companion and r/ is a parameter of order unity which 
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can be empirically adjusted to fine-tune the strength of the magnetic braking. We note 
here that the magnetic braking is effective only when the star has a considerable magnetic 
activity. It has been suggested that magnetic activity would suddenly decrease drastically 
when the star becomes fully convective. This happens when the mass of the star falls below 
~ 0.3Mq. Thus the magnetic braking is expected to turn off when the companion mass 
falls below this limit, and this is believed to be the reason for the observed period gap in 
CV systems. 

Some authors ha ye suggested lower strength s of magnetic braking than in the original 



formulation (see, e.g., 



van der Sluys et al. 



(|2005[ )). implemented by either (a) reducing the 



braking by a fixed fraction [i] = 0.25, say), or (b) using differ ent functional forms i n two 



Podsiadlowski et al. 



(12002[ l 



different regions of angular velocity. It was also suggested by 
that this reduction in strength may depend upon the mass of the convective envelope. The 
calculational scheme for collective properties of LMXBs which we will be describing in 
section [3] is capable of handling all such prescriptions of magnetic braking. However, we are 
interested here in determining the dependence of the collective properties of LMXBs on the 
very basic ingredients of the process of formation and evolution, and not on their details. 
We therefore assume the simple form given in eqn. IHin this work, postponing the study of 
more complicated effects to future works. 

A detached pre-LMXB will shrink through the loss of angular momentum which 
is caused by gravitational radiation and magnetic braking, described in Eq.([2]) by 
{J / J)nml = {j / J)gw + [j /J)mbi together with Eqs.dl]) and (jl]). We note that angular 
momentum loss due to gravitational waves as well as the magnetic braking depends on the 
instanteneous binary parameters, i.e., the two masses, the orbital separation and the radius 
of the companion for magnetic braking. Consequently, the evolution rate of the parameters 
of the system at any instant is completely determined by their current values. This fact 
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cal! be used to evolve the system in a straightforward scheme till the point of Roche lobe 



contact, which is given by Rc = ar^, w here for t 



le Ro che lobe radius we use the Eggleton 



prescription, as described in Paper II flEggleton 



19831). 



2.2. The bifurcation period 

Two main channels possible for LMXB evolution under the above conditions have been 
identified for a long time now. When the binary is sufficiently close to begin with, i.e., when 
the initial binary period is shorter than the so-called bifurcation period, the companion fills 
its Roche lobe while it is still on the main sequence, mass transfer starts, and the LMXB 
turns on. Further orbital evolution in this channel, which is followed by close LMXBs and 
CVs, is determined by angular momentum loss, and the orbital period decreases during 
this evolution. On the other hand, if the initial binary period is longer than the bifurcation 
period, the companion fills its Roche lobe after it finishes its main-sequence evolution and 
starts ascending the giant branch. In this channel, which is followed by the wider LMXBs, 
further orbital evolution is basically determined by nuclear evolution of the companion, 
and the orbital period increases during this evolution, which explains why it is called the 
bifurcation period. 

Detailed calculations of this bifurcation per iod (which depends on the stell ar masses 



Pylyser &: Savonijd (11988 



3), 



19891), showing 



and weakly on other parameters) were given by 
that this period lies in the range 14 -18 hours for the range of masses and parameters 
relevant for LMXBs. In the following subsections, we describe our calculations of LMXB 
evolution in the above two channels. 
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2.3. Evolution of LMXBs with main-sequence companions 

We consider first the evolution of LMXBs with main- sequence companions and initial 
orbital periods below the bifurcation period. Once Roche-lobe contact is established, mass 
transfer begins from the low-mass main-sequence companion to the neutron star. Mass 
transferred through the Li-point carries a specific angular momentum comparable to that 
of the binary orbit. Therefore, this mass cannot directly accrete onto the neutron star but 
rather forms an accretion disk around it, from where it is slowly transferred to the neutron 
star on a viscous timescale. We in this work do not consider the detailed dynamics of this 
accretion process, as it is not relevant for our purposes. All we need for our calculations 
here is the fact that the average accretion rate onto the neutron star (i.e., averaged over 
any fiuctuations caused by accretion-disk processes) is 7M, which determines the average 
X-ray luminosity. 

The condition for Roche- lobe contact described in the previous section, i.e., Rc = ar^, 
assumes in effect an infinitely sharp boundary of the companion star at its formal radius, 
where the stellar density suddenly falls to zero beyond this radius. Such a strict condition 
is of course unphysical and can be relaxed to accom odate a stella r atmosphere of rapidly 
falling density. We here follow the prescription by iRitterl (119881 ) of an exponentially 



decreasing mass transfer rate for larger separations, given by: 



M = Mo exp [^^^^^j (5) 

Mo here is the rate of mass transfer exactly at the start of full Roche-lobe overflow 
and H is the scale height of the atmosphere. We take H/Rc = 0.005 for main-sequence 
companions and H / R^ = 0.01 on the giant branch. Note that we set the modiflcation factor 
to unity in case of a slight overfilling of the Roche lobe, in order to avoid the very high. 
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umphysical rates of mass transfer which would formally result from the above exponential 
prescription, and which merely reflect the fact that this prescription is not valid when the 
Roche lobe is overfilled. This method of handling the beginning of Roche-lobe overflow 
smoothly is a natural way of incorporating low-luminosity systems in our scheme. We 
discuss this further in Sec J3.3[ 

After attaining full Roche lobe contact, mass transfer begins at its full rate from the 
companion to the neutron star. The response of the orbital separation to the mass transfer 
is dependent on the mass ratio q. It is a well-known fact that, if there is no mass loss (i.e., 
7 = 1), then for g < 1 the orbit expands whereas for g > 1 the orbit contracts. Thus the 
mass-transfer effects strengthen those of the angular momentum loss for g > 1, while for 
g < 1 the two effects oppose each other. 

We also note that, for sustained mass transfer, the radius of the star must always equal 
the Roche- lobe radius. This condition is automatically self-sustained for g < 1, which can 
be easily seen in the following way. For Roche-lobe underfilling systems, since no mass is 
transferred , the orbit shrinks due to angular momentum loss until Roche-lobe contact is 
re-established, while for Roche-lobe overfilling systems, heavy mass transfer occurs and 
the orbit widens until exact Roche- lobe contact is re-established. For g > 1, the situation 
is more complicated, but still viable unless g has a relatively large value (see below). For 
Roche-lobe underfilling systems, the same argument as above applies. But for Roche-lobe 
overfilling systems, mass transfer further shrinks the orbit and likely leads to runaway mass 
transfer and heavy mass loss from the companion. It is believed that this eventually brings 
the companion back into exact Roche-lobe contact. However, at sufficiently high values of g, 
even if the entire mass lost by the companion is lost from the system (7 = 0), the orbit does 
no t expand, and the above ar gument fails. This situation is similar to the one described 



by 



Podsiadlowski et al. 



(|2002[ ) for companion masses > AMq. Since such situations do 
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not appear relevant for LMXBs, however, we do not consider them any further here, and 
assume that the system is in constant Roche lobe contact. This condition is twofold. It 
implies that (a) the radius of the companion is equal to its Roche lobe radius, and (b) the 
rate of change of its radius also equals the rate of change of the Roche lobe radius. The 
latter condition can be rewritten as 



# = (6) 

Rc a tl 

For a complete description of the evolution, a mass-radius relation for the companion 
is also required. Since we assume the companion to be on the main sequence in this channel 
of LMXB evolution, we can take a simple power-law relation oc M". We note that the 
mass-radius realtion for the companion could, in principle, evolve because of the nuclear 
evolution of the low-mass companion as the evolution of a given LMXB system proceeds, 
but this effect is unimportant for the following reason. In this channel, LMXB evolution 
occurs on the timescale of angular-momentum loss from the system, which is shorter 
than the nuclear timescale of the companion. Further, as mass transfer proceeds and the 
compaion's mass decreases, its nuclear-evolution timescale becomes longer and longer, 
so that the extent of the companion's nuclear evolution over the entire LMXB evolution 
becomes quite tiny. Thus, nuclear evolution of the c ompanion is in effect "frozen", to use 



the succint description of 



Pylyser fc Savonijd (119881 ). over the whole LMXB evolution of 



interest to us here. Accordingly, we assume the above mass-radius relation to be static. 
In the computations reported in this paper, we use = 1, which is adequate for our first, 
approximate study here, which is intended to serve as a proof of principle. Other, slightly 
different, values of n are also discussed later when appropriate occasions come up. 

We emphasize here that the above assumption of a static mass-radius relation would, 
of course, be completely inadequate for giant companions. We give an account of our 



- 14 - 



handling of the time dependence of the mass-radius relation on the giant branch in Sec J2.4[ 
However, in this part of the work, this relation is static, so that the change in the radius 
of the mass-transferring companion is only through the change in its mass, and for the 
power-law mass-radius relation given above, we have (Rc/Rc) = —n{M/Mc), remembering 
the earlier relation between M and Mc- 

The rate of change of the Roche lobe radius can be computed by differentiating the 
Eggleton expression for the effective Roche lobe radius given in Paper II. We obtain 



dri rl [21n(l+p) 

'dq ~ 1.47p4 



(7) 



p I + p 

Here, p = q^^^. Equations ([6]) and (J?!) can be susbstituted in Eq.(l2]) with the above 
mass-radius relation to obtain the relation between the rate of mass transfer and the rate 
of angular momentum loss. This relation is: 

M 

NML 



(l-7)(l + 2/i) q n l + 'yqdlnrL 

7a — 1 H -; 

2 1 + g 2 2 dlnq 



Equation([H]) gives the rate of mass transfer in terms of the current system parameters. 
M and d are also related to each other by Eq.Q, with the aid of the above mass-radius 
relation. Thus the evolution of a LMXB with a main sequence companion is completely 
specified by the above equations for given values of the parameters 7 and /i. 

We take /i = 1 for main-sequence companions in this part of the work. A determination 
of 7 is not very straightforward, but we can proceed in the following way. First, we shall 
assume in this work that the value of 7 is as high as possible, subject to an absolute 
upper limit of unity. This assumption effectively means that the system's dynamics and 
hydrodynamics are such that it prefers channeling the mass through the Li point rather 
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than sending it out of the system as far as possible, consistent with all laws of motion. 
This is, of course, a simple assumption of convenience in a first study, whose effects can be 
studied in later, detailed studies. The second constraint comes from the fact that there is an 
upper limit to the accretion rate on the neutron star, which corresponds to the Eddington 
luminosity for the mass of the neutron star. Finally, Eq.(l8]) requires giq]^,^) < 0. This 
poses an additional constraint on 7, given by 

1 + 2/i g drL \ n 1 + 2/i g q dri 

2 1 + q 2rL dq J 2 2 1 + g 2ri rfg ^ ^ 

We apply the above constraints to obtain a consistent value of 7, which we use to calculate 
the evolution of the system with the aid of Eq.®. 

The essential features of the evolution of such LMXBs (and the related CV systems) 
with low-mass main-sequence companions, where angular-momentum loss from the system 
drives the mass transfer and the evolution, are well-known. Consider first the evolution 
of the LMXB orbit, as shown in FigH] for a system with an initial companion mass of 
Ms = 0.7Mq and an initial separation of aj = S.ORq, which corresponds to an initial orbital 
period of P*j,f, ~ 10 hr, i.e., shorter than the bifurcation period. As time progresses, the 
LMXB orbit shrinks and the orbital period decreases, even as the companion mass decreases 
due to mass transfer. Magnetic braking turns off when the companion becomes fully 
convective at a mass Mq ~ 0.3Mq. Beyond this point, the system continues evolution due 
to angular-momentum loss through gravitational radiation alone, until hydrogen burning 
is extinguished in the companion below a critical mass Mc ~ O.IM©. Subsequently, the 
companion follows the mass-radius relation for degenerate stars, the orbit expands after 
passing through the so-called period minimum, and the final product is a close binary 
consisting of a recycled neutron star and a low-mass He white dwarf. For our purposes here, 
it is sufficient to follow the main-sequence phase of the companion, and accordingly we 
terminate our computations when the companion mass reaches Mc = 0.1. As a consequence. 
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we do not follow the system quite upto the period minimum in FigUl but otherwise 
the period evolution is ve ry similar to that given in preyious works (see, e.g., Fig. 3 of 



Pvlvser fc Savoniiel fjlQSSh . Fig.3 of 



Podsiadlowski et al. 



( I2OO2I ). top panel, blue curve), and 



the magnetic-breaking turn-off point shows up as a "kink" in thi s diagram (at t ^ 1.5 Gyr 



in Figil]), as it does in the above diagram of 



Podsiadlowski et al. 



CO 




Fig. 1. — Angular-momentum-loss driven orbital evolution of a typical LMXB with a main- 
sequence companion (see text). 



Crucial for us in this work is the evolution of the accretion luminosity of the LMXB, 
which we display in Figj2] for the system whose orbital evolution is shown in FigJTJ 

The origin on the time axis in this figure is the formation time of the pre-LMXB, as 
before. The system comes into Roche-lobe contact at t ~ 0.4 Gyr (which shows up as 
another kink in FigH]), at which point the LMXB turns on, giving the initial luminosity 



-17- 




Fig. 2. — Angular-momentum-loss driven luminosity evolution of a typical LMXB with a 
main-sequence companion (see text and Fig\^. 

spike in Figj2l As time progresses beyond this, there is a sharp dechne in the luminosity 
until the magnetic braking turn-off point is reached at t ~ 1.5 Gyr (see above), at which 
time the accretion rate and so the luminosity drops dramatically, since the dominant 
source of angular-momentum loss is suddenly turned off. During the subsequent evolution, 
gravitational radiation is the only mechanism of angular-momentum loss, and the luminosity 
rises very slowly with time. 

The above luminosity profile (i.e., L — t relation) i s in general agreement with those 



Podsiadlowski et al. 



120021 ) work 



in previous works, an example of which is in the above 
(bottom panel of their Fig. 3, blue curve). Note that, since time is plotted on a logarithmic 
scale by these authors, the initial rise in the luminosity is actually resolved in their figure, 
while in our Figj2l where we have plotted time on a linear scale, it appears as a luminosity 
spike, as it does in other previous works with a linear time plot. The basic arguments which 
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determine the shape of this lum inosity pro f ile are interesting, and go back to two scaling 



laws originally demonstrated by 



Patterson! (119841 ) in a related but different context. We 



briefly recount these arguments here, as their importance does not seem to have been fully 
appreciated in the LMXB literature. 

The L — t profile, or equivalently the M — t profile, is computed by relating the relative 
rate of change of the companion mass, M /Mc-, to the relative rate of change of angular 
momentum, J/ J, as detailed in SecJH For clarifying the basic scalings, we shall give a 
simpler calculation here, wherein we shall use somewhat simpler prescriptions than in the 
detailed computations of Sec|2|, and obtain analytic results in the two regimes where angular 
momentum loss by either gravitational radiation (the GR regime) or magnetic braking (the 
MB regime) dominates. 

We note first that in the limit 7 = 1, i.e., no mass loss from the system (which is a 
good approximation for the angular-momentum-loss driven systems with main-sequence 
companions considered in this section, as our detailed computations have shown, as well as 
previous works), Eq.([2]) reduces to: 

NML 

In the GR regime, the left-hand side of the above equation is given by Eq.(|3]), while in the 
MB regime, it is gievn by Eq.(|l]). The second term on the right-hand side of this equation 
is basically the relative rate of change of the orbital radius, which can be related to the 
relative rate of change of the companion mass, -^7-, with the aid of two conditions. The first 
is the Roche- lobe filling condition, whose differential form is given by Eq.Q, and the second 
is the mass-radius relation, oc M", whose differential form, {Rc/Rc) = —n{M/Mc) was 
also given earlier. When this is done, we obtain 



J\ _ M 
j) ~ Me 

/ NML 



n \ f l + odlnrx 
- + 1 +<^g + 



2 J 2 d\nq 



:ii) 
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which is, of course, the simphfied version of Eq.® for 7 = 1. 

We now introduce a further simphfication which makes the calculation transparent 
and simple for our purposes of scaling demonstration in this part of the paper. Instead 
of the more complicated, but more accurate, Eggleton prescription (see Sec l2.1l and Paper 
II), which we used throughout our detailed computations described earlier, we now use the 
simpler, original Paczynski (1971) prescription, which is still a good approximation for the 
low compaion masses involved here, namely: 

1/3 

1 + g 



^L^iA-] (12) 



This simplifies the calculation, yielding 

dlnvL 1 1 



dlnq 3 1 + g 
which immediately converts Eq. (fTTl) into a very simple form 



(13) 



J\ M f5 n , 

7 + 77 - 9 h (14) 



J I Mc\6 2 

NML 



which is useful for our purposes here. If we take n = 1 for lower-main-sequence stars, the 
quantity within the parentheses on the right-hand side of the above equation reduces to 
(4/3 -g). 

We can now obtain our analytic approximations in the GR and MB regimes by 
equating the left-hand side of Eq.f fT4|) to the right-hand sides of Eqs.([3]) and (j4]) respectively, 
and re-expressing the stellar masses in terms of the total mass Mt and the mass-ratio q in 
the ensuing calculations (this is useful because Mt remains constant in these approximate 
calculations due to our assumption of 7 = 1 as given above). After some straightforward 
algebra, we get in the GR regime a profile: 

10-12n 
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where the time dependence comes from the fact that both a and q in the above equation 
are time-dependent, both decreasing with time, and the a — t profile being given in FigHJ 
In the special case n = 1 (see above), the profile is: 



-1 



" "^^^ " (l + g)(l-0.75g) - ^''^ 
The profile in the MB regime is obtained in a similar way, and found to be 

_12^ (1 _j_ „)5{n-l/3) 

Lmb oc Mmb oc a3n-i / \ y. (17) 

In the special case n = 1 (see above), the profile is: 

^(l + g)io/3 



Lmb oc Mmb oc a*" 



;i -0.75g)' 



The basic physics underlying the luminosity profile in Figj2] is clear now from the 

scalings evident in Eqs.f ll6p and f lTSj) . Both the GR and the MB profile factorize into 

a-dependent and g-dependent parts, and the profile is largely determined by the a-dependent 

part, as the g-dependent part changes relatively slowly (compared to the a-dependent part) 

as the low-mass companion loses mass and q decreases with ongoing mass transfer. Of 

course, since the a-dependence is not as strong for GR case as it is in the MB case, some 

modification does come in from the g-dependence, but the qualitative argument still holds. 

The sharp decline in L in the MB-dominated regime between the LMXB turn-on and the 

MB turn-off at t ~ 1.5 Gyr is basically due to the very strong a^-dependence of the MB 

luminosity in a shrinking orbit, as seen in Eq.( |T8l) . Similarly, the slow rise in L in the 

GR-dominated regime beyond the MB turn-off is basically due to the a~^-dependence of the 

GR luminosity in a shrinking orbit (duly modified by the g- dependence), as seen in Eq. (fT6|) . 

I 



These scalings were first pointed out by iPattersoru (1l984j ) in a study of CV and LMXB 



evolution, for the purpose of comparing the expected M — Porb correlation with the data, 
largely on CVs. (The transformation between a and Porb is, of course, readily obtained 
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through Kepler's third law.) This amounts to looking at a "snapshot" at the current epoch 
of a collection of systems with various stellar masses and at var ious stages of evolution, and 
comparing this data with the theoretical scaling. As iPattersonl noted, a reasonable account 



of the data was indeed given by the above scaling, which was subsequently confirmed 



by iPylyser &: Savonijd ( 1l988l ). by comparing the results of their detailed evolutionary 
computations with the same data. We have argued above that the same scaling must also 
be necessarily applicable to the evolution of a single system as time proceeds, and have 
demonstrated that this is indeed so. Accordingly, we shall call this scaling the Patterson 
scaling throughout this work. 

We have given above a first-principles derivation of the Patterson scaling to emphasize 
its fundame ntal and tr ansparent nature, and also because some details were different in 



the original 



Patterson 



work, e.g., this author used a slightly different magnetic-braking 
model. However, we find a close agreement between Patterson's result and ours, showing 



Patterson 



the robustness of the scaling. In particular, for the specific value of n = 0.88 that 
adopted, our scaling give Mgr oc a~°'^^ cx P^^ "^^, while Patterson gives Mgr oc P^b "^^- 
Similary, in the MB regime, our scaling gives Mmb oc oc Porl^, while Patterson gives 
Mmb oc Pg^l^. It is clear that the Patterson scaling is at the heart of angular-momentum- loss 
driven luminosity evolution of LMXBs with main-sequence companions found in this work 
and numerous previous works in the subject. 



2.4. Evolution of LMXBs with giant-branch companions 

When the pre-LMXB is not sufficiently close according to the criterion described 
earlier, i.e., the binary period is longer than the bifurcation period, the companion finishes 
its main-sequence evolution before Roche lobe contact. In such a case, the LMXB phase 
turns on when the companion goes into the giant branch, starts expanding rapidly and 
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overfills the Roche lobe. The actual transition timescale from the main sequence to the 
giant branch is that of the traversal of the Hertzsprung gap in the HR diagram, which is 
extremely short compared to the typical timescale of LMXB evolution, and hence can be 
considered practically instanteneous for our purposes here. The orbit cannot expand in step 
with this very rapid process to accommodate the expanding giant companion, so that large 
mass losses occur from the system, until the Roche-lobe filling condition is satisfied again. 
We assume in this work that there is no accretion during this short phase, i.e., we set 7 = 
in it. The angular momentum loss due to the non- mass- loss mechanisms are also negligible 
due to the short duration of this phase. Consequently, the evolution of the system can be 
described by integrating Eq.([2]) under these assumptions, which gives the condition: 



Here, the superscripts i and / denote the initial and the final state respectively. This 
equation is to be used along with the condition that the orbital separation after this phase 
is just enough to fit a giant donor within its Roche lobe, i.e., = Rgb{M[) / r ^{Ml /Mjqs)- 

In such LMXBs with giant companions, the mass-transfer rate is determined by the 
rate of expansion of the companion due to its nuclear evolution. We assume in this work 
that the condition of constant Roche-lobe filling is still valid on the giant branch and Eq.(|6]) 
holds for these systems as well. However, the assumption of a static mass-radius relation 
is not valid on the giant branch, and the time evolution of the companion radius needs to 
be determined in the following way. The radius of a star on the giant branch depends upon 
its mass as well as its luminosity. Now, the luminosity on the giant branch depends upon 
the core mass, which, in turn, depends upon the initial mass rather than the current mass 
of the companion. Thus the evolution of a LMXB with a giant donor depends upon the 
initial values of the parameters as well as their current values. With the aid of Eq.(48) of 




(19) 
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Hurley et al.l (|2000[ ). we can write the radius a star on the giant branch as: 



0.76^ 



= R(.g = l.lM-"-^(L^-^ + 0.383L;^ 



Here, is the instanteneous luminosity of the companion. If we define f{Lc) = 
(L^'"^ + 0.383L° ''®), rate of change of the stellar radius can be given as 

Rc M 1 df ■ 
— = 0.3 — + —Lc 

R, M, fdL, 

Analytic al fitting formulae for numerical stellar-evolution results provided by 



(20) 



(21) 



Hurley et al. 



can be used to calculate Lc and Lc, which, in turn, give the rate of expansion 
of the stellar radius. Eqs.dS]) and ([7]) can now be inserted, along with the new mass-radius 
relation given by Eqs. fl20|) and f l2T|) . into Eq.([2]) to obtain the equation governing the LMXB 
evolution for giant donors. We get: 



J 

^ ' NML 



7g - 1 + (^-^)(^ + ^^\ ^ + 0.15 + + 



1 + q 



2r, 



dq 



(22) 



Values of fi and 7 may be different here from those for the main-sequence donor case. 
Whereas the value of fi is completely unknown as per current underst anding of the subject 



it is customarily taken as unity in studies of populations of LMXBs (jPodsiadlowski et al. 



20021 ). We argue here, however, that the actual value in case of giant donors may well be 
less than this, since the outer envelope, which is the actual supplier of the transferred mass, 
is less tightly bound in the case of giants. We therefore adopt /x = 0.75 in our calculations. 

The value of 7 is subject to similar constraints as in the case of main-sequence donors. 
However, we note that in case of giants some mass loss would be inevitable. Therefore the 



upper limit on the value of 7 is un 



value is generally taken to be 0.5 (jPodsiadlowski et al. 



ikely to be unity in the c a se of giant companions . This 



2002 



Belczynski et al. 



20081 ) ■ but 
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this is a somewhat arbitrary assumption, as noted by iPodsiadlowski et al.l . We therefore 
keep this as a free parameter and study its effects on the evolution of LMXB systems. Of 
course, the other two constraints described earher, namely, (a) that due to the Eddington 
limit on the mass accreting onto the neutron star, and (b) that expressed by the condition 
/i(g;7,/i) < 0, still apply. The latter condition can be written explicitly as: 

7 (q - i±^^ + < 0.85 - i±^^ _ (23) 

^ V 2 1 + q 2rL dq J 2 1 + q 2rLdq ^ ^ 

The value of 7 is obtained subject to the three conditions given above, and it is then used 
to compute the evolution of the LMXB, which is described by Eq.( l22l) . We emphasize 
here that a knowledge of the current values of the system parameters is not sufficient for 
calculating LMXB evolution with giant-branch donors, which makes the computations more 
complicated. 

LMXBs with the giant donors end their evolution when they run out of the supply of 
transferrable mass. This means that the mass transfer basically continues until the entire 
envelope of the companion is transferred to the neutron star, leaving only its core. We can 
write this condition as = Mg^c where M^^c is the core mass of the companion, which 
depends upon its initial mass. This He-core now contracts basically like an isolated He-star, 
so that Roche-lobe contact is lost and the LMXB is turned off. We stop our numerical 
computations at this point. The final product of such evolution is a wide, circular binary 
consisting of a recycled neutron star and a low-mass white dwarf with Mwd < O.45M0 or 
so. 

Figure m depicts the luminosity evolution of such a prototype system, the corresponding 
orbital evolution being shown in FigJSl In this case, the initial secondary mass is Mg = 1.3 
and the initial orbital separation is = 11, corresponding to an initial orbital period of 
^orb ^ 2.6 days. In sharp contrast to the earlier angular-momentum-loss driven case, here 
the orbit expands continually as evolution proceeds, as found in previous works (see, e.g., 
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Fig. 3. — Nuclear- evolution driven orbital evolution of a typical LMXB with a giant compan- 
ion (see text). Line-style of curves coded by the values of •jmax cls indicated. 



Figs.2b a nd 2c of iTauris fc Savoniia fil999il . Fig.ll of iPodsiadlowski et all mm . Figs.9 



and 10 of 



Belczynski et al 



( 120081 ). which are Star Track studies of the above 2b and 2c 
cases in the Tauris-Savonije work). In these figures, zero on the time axis corresponds to 
the formation of the pre-LMXB, as before, but since the system remains out of Roche- lobe 
contact during the secondary's main sequence hfe, we have so shifted and rescaled the time 
axis as to show most clearly only the X-ray active part with non-zero LMXB luminosity, 
after the system reaches Roche-lobe contact at t ~ 4.55 Gyr. The evolutions of both the 
orbit and the luminosity found by us are rather similar to those found in previous works, 
when account is taken of the fact that the display methods in some of these previous wor ks 



are di fferent from ours, e.g., time plotted on a logarithmic scale in 



Podsiadlowski et al 



(I2OO2I ). and th e decreasing companion m ass instead of increasing time as the abscissa in 



some panels of 



Tauris fc Savonijd (119991 ). The major feature of the luminosity profile after 
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Fig. 4. — Nuclear- evolution driven luminosity evolution of a typical LMXB with a giant 
companion (see text and Figl^. Line-style of curves coded by the values of'jrnax o,s indicated. 

Roche-lobe contact and LMXB turn-on is a sharp initial decline, followed by a slower 
decline and a plateau phase, and a very slight rise at the longest times in some cases. 
An explanation of this luminosity profile involves details of stellar evolution on the giant 
branch, and so is ouside the scope of this paper, contrary to the basic, simple explanation in 
terms of Patterson scaling that was possible earlier for the angular-momentum-loss driven 
evolution of LMXBs with main-sequence companions. 
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3. X-ray Luminosity Function (XLF) Calculation Scheme 

3.1. The partial XLF 

The importance of the cosmic star formation history (SFH) and its effect on the 
pre-LMXB formation rate was noted in Paper II. It was explained there that the SFH over 
a few Gyrs previous to the current epoch plays a crucial role in determining the current 
population of LMXBs. This can also be seen from the two figures given in previous section 
depicting the evolution of a typical LMXB system with a main-sequence or giant-branch 
companion. Figures H] and S] show that, while typical evolutionary timescales as counted 
from the formation of the corresponding pre-LMXBs is rather similar for LMXBs with 
main-sequence and giant-branch companions, typically ~ 3 — 5 Gyrs, the timescales as 
counted from the turn-on points of the LMXBs are very different, typically ~ 3 Gyr 
for LMXBs with main-sequence companions, but typically ~ 0.2 Gyr, i.e., an order of 
magnitude shorter, for LMXBs with giant-branch companions. Therefore, if we observe two 
LMXB systems with the same luminosity in the current epoch, their turn-on epochs may 
differ by upto a few Gyrs, depending on the natures of their companions. In order to follow 
the evolution of LMXB popoulations, it is therefore essential to take in to consideration 
evolutionary changes in the SFR. This problem was first considered by 



Ghosh fc White 



( l200ll ). who showed that the peak in the number of LMXBs lags appropriately behind the 
SFR peak. They considered typical timescales of LMXB evolution that are comparable to 
the timescales obtained from our calculations. With the more detailed evolutionary scenario 
considered in our scheme, we are able to deal with the various timescales more accurately, 
considering them as functions of the initial parameters of the binary. This enables us to 
go one step beyond the calculations of these authors, and explore the evolution in the 
distributions of the essential collective properties of LMXBs. 

The evolution of a single LMXB detailed in Secj2] gives the complete evolutionary 
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history of a given LMXB with a specified initial state, which is determined by the initial 
mass of the companion (M* = Ms, see Paper II) and the initial orbital separation {a\ 
for brevity we drop the superscript from here on, and call it just a). In other words, the 
parameters of the system, e.g., luminosity, orbital period, instanteneous companion mass, 
and so on can be computed at any subsequent time. With this input, we can now calculate 
the X-ray luminosity function (XLF) of LMXBs in the following way, and we note in 
passing that the same method can, in principle, be applied to the distribution of any other 
collective property of LMXBs. 

Consider first all pre- LMXB systems formed with the same initial parameters. The 
evolution of all these systems will be similar to one another and can be computed with 
our scheme, giving the luminosity as a function of time. This relation between luminosity 
and time can then be inverted, so that, given the luminosity, we can determine the time 
required to attain it. Thus, if a given LMXB with a luminosity L is observed at a current 
time t, we know its formation time tf — t — tn,, where tn, is the look-back time, i.e., the 
time taken by the system since its formation to reach this value of luminosity. Now, the 
origin of time is set in these calculations at the formation of the pre-LMXB. Therefore, tf 
is the formation time of the corresponding pre-LMXB. Thus, the number of systems with 
a luminosity L, observed at present, is equal to the number of pre-LMXBs formed at time 
tf before, as described in Paper II. Since the formation rate can also be computed for the 
given set of initial parameters, we can calculate the number of LMXBs at each value of 
luminosity. By suitable binning, then, we can numerically obtain the XLF of these systems. 
We note here that this XLF is calculated for the systems which had a specific value of 
and a to start with. We therefore name this the partial XLF. 
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3.2. The full XLF 

The partial XLF, as described in the previous subsection, can be obtained for any 
initial- value set (Mg,a). The full, or integrated, XLF can then be obtained by integrating 
over Ms and a, or, in case of discrete data, by summing over all allowed values of Mg and a. 
We note here that the formation rate of pre-LMXBs given in Paper II already includes the 
distribution of Mg and a, so that one does not require further weighting in the process of 
such summation/integration. The range of these variables relevant for the XLF calculation 
was also discussed in Paper II. We simply state these ranges here, and refer to Paper II 
for further details. The range of is taken to be [0.1, 2.5]. We note here that this is the 
initial value of the companion mass, so that the instanteneous value of this mass (denoted 
in this work by Mc) will be lower than this at some intermediate time of evolution, after 
mass transfer starts. Systems with higher companion masses are expected to have higher 
accretion rates and lower values of 7, implying significant mass loss in some systems. When 
the systems are actually observed as LMXBs, companions more massive than IMq are 
rather unlikely to be found. Nevertheless, our scheme does allow higher-mass companions 
for completeness, so that systems hke Her X-1 can also be included in the scheme. 

The range of a is [1.0, 20.0]. Very few systems wider than this can come into Roche- lobe 
contact in a Hubble time. Many systems with companion masses at the lower end of the 
mass-range cannot come into Roche lobe contact at separations > lOi?©. We note here 
that the systems with very low-mass companions cannot reach the giant phase because 
their main sequence lifetime is longer than the Hubble time. This lower limit on the initial 
companion mass can be as high as O.9M0, depending upon the value of Hubble constant 
assumed. We note that these constraints do not need to be posed explicitly on the systems 
under consideration, as our scheme of computations is so designed as to automatically take 
care of such issues related to timescales. 
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The XLF thus obtained depends upon (a) system parameters at the pre-LMXB stage, 
which affect the resuhs through the formation rate of these systems, and (b) various 
parameters that affect the further evolution of LMXBs. A complete parameter study is 
outside the scope of this work. We therefore fix the values of nearly all the parameters 
which we list now. At the pre-LMXB stage, we take the CE-parameter a\ = 1.0, and the 
metallicity of the primary z = 0.02. To calculate the formation rate of the pre-LMXBs, we 



assume a M adatt-profile 



p = 4.6 (see 



Blain et al. 



br th e SFR, which is a peak-type profile with Zmax = 0.39 and 



(Il999[ ) and also relevant references given in Paper II for the types 
of SFR profile and the parameters used to describe them). We study the final XLF for both 
values of (3 introduced in Paper II, viz., the uniform and the falling-power-law distribution 
for the mass ratio, and also for the with kick and without kick scenarios (this latter scenario 
was shown in Paper II to give results almost identical to those for ECSN-type low kicks). 
We note here that z = 0.02 may be an incorrect assumption, since many LMXBs are found 
in early-type galaxies, where the progenitor stars could have been metal-poor. However, we 
do not test this assumption at this stage, noting that the effect of changing the metallicity 
can be judged from the variation of the pre-LMXB PDF calculated in Paper II. Values of 
the other parameters adopted in the calculation of the evolution of the LMXBs are given in 
Secs J2.3l and 12.41 We mention here again that we take fi and 7max as unity for the evolution 
with main-sequence companions, where jmax is the upper limit for 'j. fi = 0.75 is taken for 
the evolution with companions on the giant branch, and in this case we study two values of 
Imax which are 0.5 and 0.8. We realize that these particular choices allow us to explore only 
a limited section of the full parameter space. However, our aim here is to demonstrate that 
the straightforward scheme described in this work can be used to calculate the XLF, which 
can be directly compared with observations, i.e., a proof of principle. We thus concentrate 
in this work on an exploration of the relative importances of various physical processes on 
the LMXB XLF, deferring parameter studies to future works. 
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3.3. Properties of the calculated LMXB XLF 

We first note that, due to the calculational procedure adopted here, the cumulative 
XLF is easier to calculate than the differential one. Therefore, we discuss various properties 
of the calculated cumulative XLF in this subsection. We first discuss the dependence of 
the computed XLF on the pre-LMXB parameters and then that on the LMXB-evolution 
parameters. We chose two pre-LMXB parameters to illustrate the essential points. The 
first is the nature of the primordial g-distribution, and the second is the SN-kick scenario. 
We had noted earlier that ECSN-type kicks with a = 26.5 km/s give PDFs similar to the 
no-kick scenario, whereas ICCSN-type kicks with cr = 265 km/s result in a different PDF. 
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Fig. 5. — Dependece of cumulative LMXB XLF on pre-LMXB parameters. The four cases 
are coded by line-style according to the vaues of the parameters, as indicated. 

Figj5] shows the cumulative XLF for various values of pre-LMXB parameters, the four 
cases shown being the same as those considered in Paper II. (For definiteness, the XLF 
has been normalised to a total of 100 systems here.) The following general features can 
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be seen. The shape of the XLF is much more complex than that of the HMXB XLF. It 
starts with a flat region at low luminosities, which extends upto L ^ 10^^ erg/ s. We discuss 
the XLF in this low-luminosity regime later in this section. A power-law like fall is seen at 
higher luminosities, which extends upto ~ 10^^ erg/s. The XLF cuts off beyond this point, 
where the luminosity approaches the Eddington luminosity for neutron stars. We note here 
that, since the mass of the neutron star need not be exactly the same for all systems, the 
Eddington cut-off may also have a range in general. 

It can be seen that the "with-kick" scenario does not produce a very different XLF 
from the "without-kick" scenario, except for the sharper drop near the cutoff when the kicks 
are included. By contrast, the two different g-distributions produce quite different XLFs in 
the falling-power-law regime, with a steeper fall for the distribution of q with /3 = —2.7. 
The numerically computed XLFs in the luminosity range 3 x 10'^^ — 10^^ erg/s, as well as 
over other ranges, can be fitted by power-laws of the form A^(> L) oc L". The results of 
such least-squares fits are given in Table [H where the leftmost column shows the luminosity 
range considered, and the next columns give the values of n for the four different cases. 

Relative contributions of main-sequence and giant-branch donors can influence the 
shape of the overall LMXB XLF, since many characteristics of these two donor classes 
differ from each other. First, the giants come into Roche- lobe contact later than the 



Luminosity range 
(xio36 erg/s) 


13 = 
without kick 


13 = 
with kick 


13 = -2.7 
without kick 


/3 = -2.7 
with kick 


3 - 10 


-1.18 


-1.16 


-1.65 


-1.51 


10 - 100 


-1.41 


-1.44 


-1.7 


-1.71 


3- 100 


-1.31 


-1.28 


-1.59 


-1.53 



Table 1: Best- fit power-law indices over different luminosity ranges for our computed LMXB 
XLF, for different values of pre-LMXB parameters as indicated. 
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main-sequence companions. Thus the giant population corresponds to an earher population 
of primordial binaries. Since the Madau profile suggests a peak in SFR at 2; = Zmax = 0.39, 
this effect can be important. Second, systems with giant companions are typically brighter, 
but have shorter lifetimes. This makes them much smaller in number, but situated at 
the high-L end of the XLF. We study the effect of the fraction of systems with giant 
companions, which we name the giant fraction, on the XLF in the following way. We first 
evolve the LMXB systems in the way described in Sec|2l We calculate the XLF obtained 
from this evolutionary scheme which includes main sequence as well as giant donors. We 
then modify our algorithm of evolution in such a way that, when a companion completes 
its main sequence life, the system is removed from the computations. This eliminates the 
possibility of having giant donors, keeping in our computation only those systems which 
reach Roche-lobe contact during their main-sequence life, and so ensuring that all LMXBs 
in our computation have only main-sequence companions. The XLF computed in such a 
way is thus that corresponding to main-sequence donors alone. 

Figure [6] shows the XLF for only main-sequence donors, superposed on the XLF after 
the inclusion of the giants. We see that the total number of systems with giant companions 
is not more than a few percent of the total number, but the giant-companion systems all lie 
at the bright end of the XLF, making the XLF considerably flatter at this end (see Table 
1) than that with main-sequence companions alone, which has a slope of n ^ —1.94 at this 
end. We note here that, though the computation presented here to illustrate this point is for 
P = and the "with-kick" case, the qualitative result is the same for all other combinations 
of /3-value and kick-scenario. Therefore the giant fraction, although numerically small, is a 
very important factor in deciding the shape of the LMXB XLF at the bright end. 

The parameter 7 denotes the accretion rate onto the neutron star per unit rate of mass 
loss by the companion. This parameter is close to unity for main-sequence companions. 
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Fig. 6. — Relative contributions of main-sequence and giant donors in the computed XLF. 
Cases encoded by line-style as indicated. 

but its value is expected to be lower for giant companions, and this value was taken to be 
7gb = 0.8 in previous calculations. A change in the value of •jgb results in a change of 
the luminosity of those systems with giant donors which are operating at a minimum mass 
loss from the system (see Sec l2.4p . The importance of the giant fraction in determining the 
shape of the XLF has already been shown above, from which we expect that lowering the 
value of 'Jgb will lower the influence of the giant-companion systems on the total XLF (since 
the maximum allowed luminosity for these objects will be lowered), and so make the XLF 
fall more steeply at the bright end. Figure [7| shows that this is indeed the case, displaying 
the effect for two values of 7g_b- The first value is 0.8, which has been used throughout 
this work. The second value is 0.5, which has been suggested by some authors, as noted in 
Sec J2.4[ The power-law index obtained by least-squares fit method in the luminosity range 
LsQ = 3 — 100 is n = —1.48 in the second case, i.e., steeper than in the first case (see Table 
[1]), but not as steep as in the 'Jgb = 0.8 case (see above). 
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Fig. 7. — Dependence of LMXB XLF on ■Jgb- Curves labeled by the respective values of 

iGB- 

Finally, we consider the nature of our computed XLF at very low luminosities [L < 10^^ 
erg/s). The XLF at these luminosities appears to be not very interesting for several 
reasons. First, it is difficult to observe such faint systems with the current sensitivity of 
the X-ray detectors, even in nearby galaxies. Thus the observed XLF will be strongly 
influenced by selection effects in this region. Although studies often attempt to account for 
incompleteness, there are necessarily large error-bars in observed XLFs in this regime, as a 
study of relevant works shows. Second, one can see from our calculations that the standard 
theory predicts that there would be very few systems in this range. In fact, if we exclude 
Roche-lobe underfilling companions, the minimum luminosity possible is ~ 10^^ erg/s, 
which corresponds to a binary of a IAMq neutron star and a main-sequence companion 
with a mass just below O.SMq, the critical point where the magentic braking is turned 
off (corresponding roughly to the minimum luminosity in Figj2]). Only with our inclusion 
of Roche-lobe underfilling companions and the Ritter recipe for atmospheric Roche-lobe 
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overflow (wherein the accretion rate drops exponentially for increasing disparity between 
the stellar radius and the Roche-lobe radius; see Eq.([5])), as detailed in Sec J2.3l do we 
account for systems at lower luminosity. 
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99.85 
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^ 99.75 

-z. 

99.7 
99.65 

99.6 
99.55 

99 5 

0.001 0.01 0.1 
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Fig. 8. — LMXB XLF at very low luminosities. 

Figure [H] shows the details of the computed XLF in this luminosity regime, using the 
above Ritter recipe. Because of the very narrow range of the ordinate (A^(> L)) involved 
in this regime, this part of the XLF appeared almost flat (i.e., constant A^(> L)) in all 
previous XLFs shown in this paper, where the scale of A^(> L) was logarithmic. In this 
figure, we have displayed A^(> L) on a linear scale, and shown only that range of A^(> L) 
over which it varies in this regime, in order to show the nature of its variation. Apart from 
a artifacts resulting from binning and numerical effects, the average trend is consistent with 
a logarithmic decrease in A^(> L) with increasing L-values. 
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4. Discussion 



4.1. Comparison with observations 

In recent years, LMXB XLFs have been obtained for many early-type galaxies 
from Chandra and XMM-Newton observations. Many such works have come to the 



conclusion that the observed LMXB X LF cannot be described by a simple, sing 



distribution, unlike th e HMXB XLF ( 



Kim fc Fabbiano 



Fragos et al. 



2004 



Grimm et al. 



Revnivtsev et al. 



2008 



2002 



Hurnphrey fc Buote 



Voss fc Gilfanov 



2006 



e-pow er-law 



2004 



Voss et al 



12009 



20091 ). These works have argued that the LMXB XLF can be adequately 
described by a broken power-law with a high luminosity cut-off. It was suggested by 



Grimm et al 



(|2002[ ) that the number of LMXBs in a galaxy would scale with the total 
stellar mass of that galaxy, so that, when the LMXB XLF is normalised by this stellar 
mass, a roughly universal LMXF XLF would emerge, similar to the situation for the HMXB 
XLF and its scaling with the current SFR in the particular galaxy, which was described in 
Paper I. However, this scaling was not found to be as precise as that in the case of HMXBs. 



It was shown by I Gilfanov! (|2004| ) that a doubly-broken power-law (i.e., a power law 
with two breaks) template gave an adequate fit to the observed universal LMXB XLF. The 
break points obtained by these authors were at L^q = 19 and L^q = 500, and the power-law 
exponents given by them for the differential XLF dN/dL in the three regions separated 
by these breaks were —1.0 ± 0.13, —1.86 ± 0.12 and —4.8 ± 1.1, going from low to high 
luminosities. We stress again that these authors provided fits to the differential XLF, i.e., 
dN/dL. Power-law indices in this description are obtained by subtracting 1 from the ones 
in the cumulative description, except in the special case of a power-law index of —1 in the 
differential description, where in the corresponding cumulative description A^(> L) is a 
logarithmically decreasing function of L. 
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However, iKim fc Fabbiand (12004 ) proposed a different template. These authors showed 
that the XLF of a single galaxy can be described adequately by a single power-law. With a 
sample of 14 E and SO galaxies, these authors concluded that power-law exponents in the 
range (—0.8, —1.2) could explain the cumulative XLFs of all galaxies in their sample, and 
proposed further that a single power-law with an exponent —1.1 ±0.1 for the cumulative 
XLF would be consistent with the combined data. However, these authors also pointed 
out that a broken power-law would provide an improved fit . The bre ak-point suggested by 



these aut hors roughly agre e s wit 
given by 



Gilfanovl . The power-law slopes 



1 the second break-point of 
Kim fc Fabbiand (j200J) for the differential XLF were -0.8 ± 0.2 and -1.8 ± 0.6 
respectively, below and above this break. Thus, a broken power-law (with one or two 
breaks) generally seems to account for all observations of LMXB XLFs to date, and we 
shall use this as the template for our discussion here. 

However, comparing these observed XLFs with our computed XLF is not completely 
straightforward, due to the complex nature of the XLFs. Note first that observed XLFs 
include NS as well BH binaries, and so extend upto luminosities ~ 10^^ erg/s. Systems 
beyond the second break in the above template are clearly BH systems, whereas their 
contributions below the second break, which need not be zero, are unknown. By contrast, 
the calculations presented by us in Paper II and in this work include only NS systems, 
so that a strong cut-off is expected and seen at the neutron-star Eddington luminosity. 
The power-law regime in our computed XLF is in the range L36 = (3, 100), which roughly 
overlaps with the middle region between the two breaks in the observed "universal" XLF, 
which is L36 = (20, 500). We note from Table [T| that our /3 = case gives a power-law slope 



relatively close to that of the observed XLF of 



Gilfanovl (l2004j ). whereas the {3 = —2.7 case 



gives an XLF considerably steeper than t 
also match reasonably with the results of 
fit to individual galaxies. 



le observed one. Our re sults for the (3 = case 



Kim &: Fabbiand (l2004l ) for the single power-law 
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On the whole, it appears at this stage of comparison between observed and computed 
XLFs that [3 = provides a better match than /3 = —2.7. However, it remains to be 
explored if different appropriate choices of other parameters might not give a power-law 
slope in better agreement with observations even for the latter value of (3. 

4.2. Conclusions 

We remind ourselves that, while the computed XLF presented in this paper depends 
closely on the computed pre-LMXB PDFs presented in Paper II, pre-LMXBs are, by and 
large, not observable. Accordingly, the LMXB XLF comparisons presented in this paper 
(and possibly such comparisons of calculable distributions of other observable collective 
properties of LMXBs) serve as a major check on the above pre-LMXB PDFs. To the extent 
that even a qualitative account of observed LMXB XLFs is possible by our simple scheme, 
a certain measure of confidence in the inferred pre-LMXB PDFs of Paper II is gained from 
these investigations, which has implications for studies of related types of X-ray binaries. 

We saw in the previous section that the giant fraction is an important parameter in 
determining the shape of the XLF at the bright end. We emphasize that this fraction is 
not a free parameter in our calculational scheme, set by hand. Rather, it is decided by 
the shape of the pre-LMXB PDF. Thus, a pre-LMXB PDF which gives a larger giant 
fraction produces a shallower XLF. At one remove, this corresponds to the inclusion of 
a higher fraction of wider pre-LMXBs with higher-mass companions. The choices made 
by us for primordial-binary parameters and subsequent evloutionary parameters in Paper 
II determined this fraction entirely, leaving no free parameters to adjust at the stage of 
computing the XLF. 

A primordial-binary parameter which seems to be very important for its ultimate 



-40 - 



influence on the XLF shape is that which controls the shape of the ^-distribution, viz., the 
parameter (3. The reason for this can be traced to the influence of P on the pre-LMXB PDF 
for both Ms and a, as shown in Paper II (see Figs. 6 and 11 of that paper). Between the 
gr-distribution and the SN-kick scenario, the former has by far the dominant influence on the 
pre-LMXB PDF, when other parameters are held fixed. It comes as no surprise, therefore, 
that (3 should have a strong influence on the XLF shape. What wc do find interesting 
is that, of the two values of (3 used in previous works on the subject, one fares so much 
better than the other by the XLF criterion. However, we must remind ourselves that this 
has been shown to be true in this work only for the chosen values of the other parameters. 
It would be premature to take the results of this work as a definite indication that a fiat 
g-distribution in primordial-binary distribution is always preferred, until the phase-space of 
the other parameters is thoroughly studied. However, to the extent that the values of the 
latter parameters are representative of actual pre-LMXBs, our result does seem suggestive. 

By contrast, the SN-kick scenario has little infiuence on the XLF, as expected form the 
above arguments, since it has little effect on the pre-LMXB PDF, as shown in Paper II. 
Indeed, the only case in which it has any significant effect is that of the a-distribution in the 
^ = case (see Fig.ll of Paper II). But even this has little final infiuence on the XLF, since 
the Mg-distribution is much more effective in controlling the XLF than the a-distribution. 

The calculations of various rates of the mass transfer and the orbital shrinking with 
our prescription of the mass loss depend crucially upon the assumption of constant Roche 
lobe contact during mass transfer and loss (if any). We have given arguments in favor 
of this assumption at appropriate places in this work. However, we stress here that this 
assumption should be made with appropriate timescales in mind, as we now explain. The 
typical timescales of the LMXB evolution are ~ 2 — 3 Gyrs. A timestep of 10~^ Gyrs is 
therefore sufficiently short for keeping an accurate track of the evolution of these systems. 
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and we have, of course, verified that finer time-steps give essentially the same result. This 
is tantamount to assuming that the orbit readjusts itself to Roche- lobe contact within the 
timescale of the step size assumed, i.e., ~ 10^ yr, which appears to be a safe assumption. 

With the cautionary remark that, since the falling-power-law region in our calculated 
XLF does not exactly overlap with the power-law region suggested by empirical fits to 
the observations, we should be careful about drawing conclusions, it does seem significant 
that a uniform primordial g-distribution does consistently lead to a better agreement 
with observations than the falling-power-law g-distribution with a slope of -2.7. With due 
caution, we therefore suggest that our results do indicate that a /3 = distribution for 
LMXB-progenitor primordial binaries seems to be favored by the XLF observations. It 
would be interesting to see if the inclusion of black-hole sytems leads to an exact match 
with the observations, which would be a rather gratifying confirmation of our suggestion for 
the underlying primordial g-distribution. We note, however, that the other parameters held 
fixed in this study, e.g., the CE parameter, may also have significant roles in determining 
the slope in this region. The value of 'jmax on the giant branch is also uncertain. However, 
low values of this parameter around 0.5 would decrease the effect of the giant populations 
and hence are less likely. It may also be not impossible that this parameter gradually 
decreases with time as the companion star evolves along the giant branch. As we saw 
earlier, the value 'jmax = 0.8 produces a better fit with the data. This may thus represent 
an average value of the 'jmax over this entire evolution. 

The problems of a proper representation of the LMXB XLF at the lowest-luminosity 
end are many, not the least of which is the possibility that the accretion paradigm in 



this region ni ay be entirely different, e.g., ADAF (see iNarayan fc McClintockl (120081 ): 



Lasotal (l2008l ) for recent reviews) or related flow models. Incorporation of such models 



into our scheme is a very ambitious task, since the relations between the mass-transfer 
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rate, the mass- accretion rate, and the output X-ray luminosity are thought to be very 
different in some of these models from the simple accretion-disk paradigm we have used 
here. Accordingly, we defer the inclusions of both black-hole systems and such acretion 
flow models to future works. We conclude by re-emphasizing that our effort in this series 
of papers must be regarded as an attempt at a proof-of-principles exploration of whether 
observed collective properties of X-ray binaries can be accounted for by evolving canonical 
collective properties of primordial bianries through well-accepted scenarios of individual 
binary evolution. Considering the simplicity of our scheme (for both HMXBs and LMXBs), 
the agreement for HMXBs was remarkable, and that for LMXBs, although not as precise, is 
still qualitatively in the correct direction. Thus encouraged, we feel justified in attempting 
more elaborate future explorations, which would incorporate detailed stellar-evolution 
models. 
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